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We study the asymptotic long-range behavior of the time-dependent correlation function of the surface charge 
density induced on the interface between two media of distinct dielectric functions which are in thermal equi- 
librium with one another as well as with the radiated electromagnetic field. We start with a short review which 
summarizes the results obtained by using quantum and classical descriptions of media, in both non-retarded and 
retarded regimes of particle interactions. The classical static result for the flat interface is rederived by using 
a combination of the microscopic linear response theory and the macroscopic method of electrostatic image 
charges. The method is then applied to the case of a spherically shaped interface between media. 
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1 A short review 

This article is about a simple inhomogeneous physical system consisting of two distinct electrical media in contact 
with one another which are in thermal equilibrium. Such system is easily accessible in experiments and can serve 
as a tool for testing quantum mechanics (and its classical limit) and crossover between retardation and non- 
retarded effects. 

The special geometry of two semi-infinite media with a flat interface, formulated in the three-dimensional 
Cartesian space of points r = (x, y, z), is pictured in Fig. Q] The two media with the distinct frequency-dependent 
dielectric functions ei{u>) and e2(u>) occupy the complementary half-spaces x > and x < 0, respectively. The 
flat interface between media is the plane x = 0. We recall that, in Gauss units, e{uS) = 1 for vacuum, while the 
static dielectric constant e(0) is finite (> 1) for dielectrics and diverging, e(0) — > ioo, for conductors. The system 
is translationally invariant along each plane formed by the two coordinates R = (y, z) perpendicular to x. We 
assume that the media have no magnetic structure and the magnetic permeabilities [i\ = /12 = 1. Non-relativistic 
charged particles forming the two media and the radiated electromagnetic (EM) field are in thermal equilibrium 
at some temperature T, or the inverse temperature (3 = 1/(JsbT) with fee being the Boltzmann constant. 
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x 

Fig. 1 The geometry of two semi-infinite electrical me- Fig. 2 The geometry of two electrical media with the 
dia with the flat interface contact at x = 0. spherical interface of radius a. 

The different electric properties of the media give rise to a microscopic surface charge density (operator in 
quantum mechanics) er(£,R) at time t and point r = (0,R) on the interface. It is understood as being the 
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microscopic volume charge density integrated along the x-axis on some microscopic depth from the interface. 
According to the elementary electrodynamics, the surface charge density is associated with the discontinuity of 
the normal x component of the microscopic electric field at the interface: 

a(t,R) = i- [E+(t, R) - E~(t,K)] , (1) 

where the superscript +/— means approaching the interface through the limit x —> + /0~. The tangential y and 
z components of the electric field are continuous at the interface. The correlation of fluctuations of the surface 
charge density at two points on the interface, with times different by t and distances different by R = |R|, is 
described by the (symmetrized) correlation function 

S(t, R) = \ (<r(t, R)cr(0, 0) + a{0, 0)<r(t, R)) T , (2) 

where (• • • ) T represents a truncated equilibrium average at the inverse temperature /3, (AB) T = (AB) — (A)(B). 
With regard to ((TJ, this correlation function is related to fluctuations of the electric field on the interface. 

Although the system is not in a critical state, the combination of the spatial inhomogeneity and long-ranged 
EM interactions in the media causes that the asymptotic large-distance behavior of the surface charge correlation 
function (0 exhibits, in general, the long-range tail of type 

pS(t,R)~^, i?^oo (3) 

with some prefactor (slope) function h(t). It is useful to introduce the two-dimensional Fourier transform 
S(t,q) = J d 2 i?exp(— iq • R)S*(t,R) , where the wave vector q = (q y ,q z ). In the sense of distributions, 
the Fourier transform of 1/i? 3 is — 2iiq. Consequently, (3S(t, q) has a kink singularity at q = 0: /3S(t, q) ~ 
—2irh(t)q for q — > 0. S(t, q) is measurable by scattering experiments. The fact that it is linear in q for small q 
makes this quantity very different from the usual bulk structure functions with short-range (usually exponential) 
decays which are proportional to q 2 in the limit q — > 0. 

The form of the prefactor function h(t) depends on "the physical model" used. Defining media, one can apply: 

• Either classical mechanics, h c \{t), or quantum mechanics, h qu (t). According to the correspondence princi- 
ple, the quantum mechanics reduces to the classical one in the high-temperature limit. 

• Either non-retarded description, h^ nT \ or retarded description, h^'\ of particle interactions. In the non- 
retarded regime, the speed of light is taken to be infinite, c = oo, ignoring in this way magnetic forces. 
In the retarded regime, c is assumed finite and so the charged particles are fully coupled to both electric 
and magnetic parts of the radiated EM field. According to the Bohr- Van Leeuwen theorem [T|, magnetic 
degrees of freedom can be effectively eliminated from statistical averages of classical systems for the static 
case t = 0, so that h^\o) = h^ r) (0) = h cl (0). If the EM field is fhermalized, it was shown in Ref. O 
that the decoupling arises only at distances larger than the thermal photon wavelength f3Kc, where the EM 
field itself can be treated classically. Quantum electrodynamics must be used at shorter distances, which 
invalidates the decoupling. At room temperature, the wavelength (3hc is rather large compared to usual 
microscopic scales and might have some experimental importance. The decoupling is preserved in our 
macroscopic electrodynamics. For time differences t ^ 0, magnetic interactions do affect classical charged 
systems in thermal equilibrium and so, in general, h^\t) ^ n£ (t). 

The media configuration studied in the past was exclusively a conductor in contact with vacuum. When the 
conductor is modelled by a classical Coulomb fluid, a microscopic analysis J5] leads to the universal static result 

h c \(0) = (conductor vs. vacuum), (4) 

independent of the fluid composition. The same result has been obtained later |4'| by using simple macroscopic 
arguments. As concerns quantum microscopic models of Coulomb fluids, a substantial simplification arises in 
the so-called jellium (one-component plasma), i.e. the system of identical pointlike particles of charge e, mass 
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m and bulk number density n, immersed in a fixed uniform neutralizing background of charge density —en. The 
dynamical properties of the jellium have a special feature: There is no viscous damping of the long-wavelength 
plasma oscillations [5]. In the non-retarded regime of the Maxwell equations, the frequencies of nondispersive 
long-wavelength collective modes, namely lu p of the bulk plasmons and lj s of the surface plasmons, are given 
by uj p = \J 'Anne 2 / m and u> s = ui p / \/2. The dielectric function of the jellium is described by a one-resonance 
Drude formula 

e(w) = 1 - 



uj{u) + irj) ' 



(5) 



where the dissipation constant r\ is taken as positive infinitesimal, r\ 
has the non-universal form |6| 



+ . The obtained time-dependent result 
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/3h J ' 



(6) 



(7) 



In the complex upper half -plane, g(u>) possesses an infinite sequence of simple poles at the (imaginary) Matsubara 
frequencies i£j (j = 1,2,...). In the classical limit f3focj — > 0, the function g(ui) — > 1 for any u) and the quantum 
formula ((6} reduces to 



1 

8^ 



2 L 2 cos(uj s t) — cos(wpt)] . 



(8) 



For t = 0, we recover the universal classical static result (|4). 

To study a general configuration of two media in contact and the effect of retardation, in a series of works []7]|8] 
we applied the macroscopic fluctuational electrodynamics of Rytov (9) to the present problem. The classical static 
result was found in the form 
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where t\ = t\ (0) and = £2(0). For the previously studied configuration of a conductor (ei — 
(e2 = 1), we recover the universal formula (3). In the retarded regime, for any t we found that 

/&>(*) = Mo), 



(9) 

ioo and vacuum 
(10) 



independently of %, c and the temperature. This surprising result was verified explicitly on a microscopic model 
of jellium with retarded interactions fl8]. In the non-retarded regime, we got 



/i (nr) (i) = — 
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where f(it) is an analytic r 

f(r) = 2 
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cos (^r) 



(11) 



(12) 



Plugging into this expression the dielectric function of the jellium (0, one recovers the microscopic result (|6). A 
crossover between the retarded and non-retarded regions was found at q cl - oss ~ w p /c: If g < Across (<? > 9cross)> 
the slope function ft,( r ' (/i^ 11 ^) takes place. 

After the short overview, we shall investigate how the classical static result <j9j modifies itself for a spherically 
curved interface between two media (see Fig. [2}- Firstly, in section 2, we provide an alternative derivation of ((9) 
for the flat interface by using a combination of microscopic and macroscopic approaches. This method allows us 
to generalize the classical static result to the spherical interface between media in section 3. 
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2 Rederivation of the classical static result © for the flat interface 

Our strategy is to compute, for an arbitrary configuration of two points r and r' in the media 1 and 2, the 
correlation function (</>(r)0(r')), where (f>(r) is the microscopic electric potential created by constituents of the 
media at point r. It is related to the microscopic charge density p by <f>(r) = J dr"/?(r")/|r — r"|, where the 
integral is over the whole space. Since E(r) = — V^(r), we have 

{E x (r)E x (r')) T = ^ 7 (0(r)0(r')) T (13) 
According to the relation ([1), the surface charge correlation function is given by 

(er(R)(7(R')) T = JL (E+ (R)E+ (R) + E~ {R)E~ (R) - E+ (R)B" (R) - E~ (R)E+ (R)) T . 

(14) 

Let us introduce an infinitesimal test charge q at point r. Microscopically, denoting by </>tot(r') the total (i.e. 
direct plus due to the particles) microscopic potential induced at point r', it holds (^tot(r'))<? = 9/l r — r 'l + 
(<p(r')) q , where (• • • ) q means an equilibrium average in the presence of charge q. The additional Hamiltonian is 
SH = qc/)(r). Using the classical linear response theory for (</>(r')) g , we obtain 

(4>{v')) q = (</>(/)) - f3(^(r')5H) T = (0(r')> - /3g(0(r')^(r)) T , (15) 
where (• ■ ■ ) means the standard equilibrium average in the absence of the test charge q. Consequently, 

/3(0(r)0(r')) T = - - q [(^ot(r')> 9 - W))] ■ (16) 

The shift of the mean potential at r' due to the test charge q at r, (</>tot(r'))<j — (</>(r')), is determined by using 
the macroscopic method of images 110111 II . If both points r and r' are inside medium 1, the shift is given by 

<**C»« - ^ - 7^7\ + ^ - (17) 

where r* = (— x, R) is the position of the image charge q' . If both points r and r' are inside medium 2, the 
formula ( fTTI i with interchanged indices 1 and 2 applies. If r is in medium 1 and r' is in medium 2, or vice versa, 
it holds 

(0 tot (r')>, - <#r')) = -r^— p q" = ~^q- (18) 
£2|r'-r| £l+£2 

Inserting the shift of the mean potential to ( [ToT i for all four possible configurations of points r and r' in media 
1 and 2, using ( fT3l and the relations 



d 2 1 



dxdx 1 |r' — r| 

we find that 



c) 2 



x—x'—O 



|R-R'| 3 ' dxdx' |r' -r*| 



x—x'—O 



|R-R'|3' 



(19) 



P(Et(R)Et(K> r = (l - J- + -|_) (20) 
Considering these correlators in fH4] . we end up with the classical static result ((9J- 
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3 Classical static result for a spherical interface between media 

The same method is now applied to the spherical interface of radius a between media 1 (r > a) and 2 (r < a), 
see Fig. [2] The microscopic formula ([Tot remains unchanged. The macroscopic shift of the mean potential at r' 
due to the test charge q at r for the spherical geometry can be expanded in terms of the Legendre polynomials 
{P;(cos 6)}^ , with 9 being the angle between r and r', as ifTTI 



oo 

- [(0 tot (r')) 9 - <«Kr')>] = $>(r,r')fi(co S 0). 



(22) 



Denoting by 7-< (r>) the smaller (larger) of r = |r| and r' = \r'\, one has 



r') 

9i(r, r') 
9i(r, r') 
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J 



(e 1 -e 2 )l a 



21+1 



< ei(l + l) + e 2 l r 1 * 1 
1 , (e a -ei)(Z + l) r 



i+i 



> 



e 1 (l + l) + e 2 l a 2 '+! 



for r > a and r' > a, 



for r < a and r' < a, 



21 + 1 



for r > a and r' < a. 



ei{l + 1) + e 2 l H+ 1 

We define the image position of a charge at r as r* = (a/r) 2 r and consider the summation formulas |[T2l 

J. 



1 



|r-r'| y^iTp 



cos( 



r' 2 - 2rr'cos6 ^ ' ' 



(23) 
(24) 
(25) 

(26) 



a 1_ 1 J Ei^o ^Ucosfl) ^'yli if r > a and r' > a, 

r\r*-r'\ ~ + {rr') 2 /a 2 - 2rr' costf ~ \ J2Zo p i( cosd )^0 if r < a and / < a, 

(27) 



VPi(cos6») =/ rW for \t\ < 1 and a > 0, (28) 

^ l + a J v /l + (te) 2 -2tecos0 

to obtain, using (fTol l. the two-point electric potential correlation functions: 

W(rW)} T = (l 



ex J |r - r'| ei(e x + e 2 ) r |r* - r'| 

+ 7^-^- /'dx - (29) 
(£1 + 62) Jo y/ a 2 x 2 + (rr 1 ) 2 / a 2 — 2rr'x cos 9 



for r > a and 7-' > a, 



/3 Wl W)} T = (1 



e 2 J |r- r'| e 2 (ei + e 2 ) r |r* - r'| 

+ 7^^ dx - (30) 

(ei +£2) Jo v fl2 + (rr'x) 2 /a 2 - 2rr'xcos9 



for r < a and r' < a, 



/3(^(r)0(r')) T = ( 1 — 1 r^—r + , 61 £2 , 9 / drc - X I (31) 

£l +£2; |r-r'| (ei+e 2 ) 2 7 \A 2 + (r'z) 2 - 2rr'x cos 
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for r > a and r 1 < a. Here, the exponent v = £2/(61 + £2) lies in the interval [0, 1]. 

The radial component of the electric field, normal to the surface of the sphere, is given by E n (r) — —d r (f>(r). 
Thus the counterpart of (O reads (E n (r)E n (r')) T = <9 2 r , (4>(r)4>(r')) T . Using the relation (fT4l) . with E± 
substituted by after lengthy algebra we find the surface charge correlation function between two points R 
and R' on the sphere interface of the form 



/?<<r(R)a(R')> T 
where the integral 



1 

"8^2 



1 

(-2 



ei - £2 1 



£1 + £ 2 J |R - R'| 3 16^ 2 (£1 + £ 2 ) 2 <z 3 



I(y, cos 6') = / dxx 



3(1 -a; 2 ) 2 



2(1 



(1 + x 2 - 2xcos8) 5 / 2 (1 + x 2 - 2xcos9) 3 / 2 



I {v, cos 6), (32) 



(33) 



Here, the angle 9 between the points R and R' on the sphere of radius a is given by cos 9 = 1 — | (| R — R'|/a) 2 . 
The above two formulas hold for all macroscopic distances |R — R'| and sphere radiuses a. 

A special case of physical interest is the fixed distance |R — R'| and the large radius limit |R — R'|/a — > 0. 



Writing 1 



2x cos 9 = (1 — x) + ex with e = (|R — R'|/a) <§; 1 and using the small-e expansions 



dx 
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(fj, > —1), we obtain cos 9) 
/3<<7(R)<r(R')) T 



_ 1 
ex] 3 / 2 e 

= 2/e + 0(1). Thus from (O we conclude that 
1/1 1 4 \ 1 

\7i + 7 2 ~ £i + £ 2 J |R-R| 3 
1 £1 - £2 1 



(34) 



(35) 



e 2 ) 



|R- R 



(36) 



Stt 2 ( £l 

The leading term, invariant with respect to the media interchange £1 o £2, corresponds to the classical static 
result for the flat interface ((9). The first curvature correction changes the sign under £1 <H- £ 2 ; this is expected 
since the media interchange is effectively equivalent to the change of the interface curvature to the opposite one. 

Another interesting case is the antipode configuration with |R — R'| = 2a, for which curvature effects are not 
erased. Substituting cos 9 = — 1 in the integral (l33l l. I(v, —1) is expressible in terms of harmonic numbers lfi"2l . 

An open question is a sphere extension of the quantum time-dependent result fllOl . derived for a flat interface 
by taking into account the retardation. This result should be recovered in the large radius limit |R — R'| <C a, 
while curvature effects might invalidate it for |R — R'| ~ a. The answer requires the application of Rytov's EM 
fluctuational theory @ to curved interfaces between media which is a difficult task left for near future. 
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